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SECTION  I 


REPORT  SUMMARY 


High  electric  fields  are  obviously  advantageous  in  powerful 
electron  accelerators:  for  a  given  energy  gain  per  particle,  the  device 
length  is  inversely  proportional  to  the  maximum  field;  and  the  maximum 
current  is  roughly  proportional  to  the  field.  The  transverse  RF  accelerator 
provides  a  high  field  strength  by  using  an  accelerator  cavity  analogous 
to  an  optical  resonator:  the  focussing  of  electromagnetic  energy  produces 
large  fields  along  the  cavity  axis  for  manageable  fields  at  the  wall.  In 
addition,  it  is  possible  to  have  the  fields  largely  transverse,  permitting 
very  high  accelerating  gradients.  This  report  discusses  two  basic  topics: 
the  field  structure  as  related  to  the  acceleration  mechanism,  and  limita¬ 
tions  on  the  electric  field  strength^.  Important  aspects  that  are  not  part 
of  the  research  objectives  are  the  accelerator  power  requirements,  the  effects 
of  cavity  loading,  and  other  phenomena. 

The  accelerator  cavity  can  be  conceptually  generated  from  a 
linear  optical  resonator  as  shown  in  Figure  la.  The  resonator  is  then 
rotated  about  an  axis  in  the  focal  plane,  to  give  the  cylindrical  resonator 
of  Figure  lb.  This  configuration  gives  large  intensities  in  the  shaded 
region,  from  the  combined  effect  of  axial  focussing  already  present  in 
Figure  la,  and  the  radial  focussing  introduced  by  the  rotation.  Finally, 
the  accelerator  cavity,  in  Figure  lc,  is  extended  axially  to  allow  wave¬ 
lengths  comparable  to  cavity  dimensions. 


i 


The  electromagnetic  cavity  modes  suitable  for  electron  acceleration 
have  a  purely  azimuthal  magnetic  field,  and  an  mostly  axial  electric  field. 

The  modes  can  be  found  exactly.  However,  in  the  interesting 
parameter  region  the  modes  are  close  to  the  well-known  higher  order 
Hermite-Gaussian  optical  resonator  modes  shown  in  Figure  2a.  The  electric 
field  direction  during  one  half-cycle  is  as  given  in  Figure  2b,  and  is 
reversed  in  the  other  half-cycle.  The  field  is  obviously  cylindrically 
symmetric,  and  extends  out  to  the  cavity  wall. 

The  accelerator  process  is  illustrated  in  Figure  3.  As  an 
electron  enters  the  first  lobe  in  the  spatial  field  pattern.  Figure  3a  , 
the  field  points  opposite  to  the  electron  velocity,  and  the  negatively 
charged  electron  increases  its  momentum  (but  not  its  velocity)  during  the 
first  half-cycle  of  the  field  o<ut<7T  .  At  wt  =  it  the  field 
reverses,  and  the  situation  is  as  given  in  Figure  3b:  at  this  time  the 
electron  shown  has  not  yet  left  the  first  spatial  lobe,  and  it  loses  momentum. 
Again,  the  electron’s  velocity  remains  constant,  so  that  a  short  time 
e/w  later  the  electron  enters  the  middle  lobe,  (Figure  3c),  where 
it  is  again  accelerated  until  it  leaves  this  lobe  or  cot  reaches  2tt  , 
whichever  is  first.  At  this  point  there  might  be  another  momentum  decrease, 
but  it  is  obvious  that  the  net  effect  of  the  various  accelerations  and 
decelerations  is  a  net  energy  gain  AE  of  order  AE  =  eEp  L/2  , 

where  Ep  is  a  typical  field  strength,  and  L  is  the  width  of  the 
spatial  pattern  as  indicated  in  Figure  3c.  From  the  discussion  it  is 
clear  that  the  acceleration  process  is,  basically,  resonant  wave-particle 
interaction  with  an  approximately  sinusoidal  wave.  The  calculated  energy 
transfer  is  shown  in  Figure  4. 

The  accelerating  field  strength  is  intrinsically  limited  by  the 
allowable  field  at  the  cavity  wall.  The  normal  field  component  will  be 
limited  by  electric  breakdown  to  30  MeV/m  or  so,  but  the  parallel 
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component  should  be  up  to  the  skin  current  heating  limit  of  =  109  V/m  (the 
estimate  is  justified  in  Appendix  I)*  Therefore  it  is  necessary  to  know 
the  direction  of  the  field  at  the  wall,  and  to  select  a  cavity  geometry  or 
accelerator  operating  regime  where  the  electric  field  is  mainly  parallel  to 
the  wall.  A  simple  argument  for  the  field  direction  gives  Ex/E (|  =  X/A.^  , 

where  A  =  u/c  is  the  free  space  wavelength  and  il||  is  a  typical  length 
scale  of  the  field  pattern  on  the  wall  (see  Appendix  II). 


The  above  estimate  indicates  that  the  large  wall  fields  predicted 
for  the  transverse  RF  accelerator  are  obtainable  when  the  free-space  wave¬ 
length  is  small  compared  to  the  wavelength  of  the  fields  on  the  wall.  This 
is  assured  if  the  cavity  is  used  with  a  high  mode  number  in  the  radial 
direction,  or  a  high  frequency.  An  approximate  formula  is  given  in 
Appendix  III. 

The  transverse  accelerator  research  reported  on  here  concentrated 
on  these  topics.  To  summarize: 

i)  use  of  laser-like  higher  order  cavity  modes 

in  approximate  resonance  with  relativistic  electrons, 
and 

ii)  the  ability  to  make  the  electric  fields  at  the  wall 
largely  transverse,  leading  to  very  high  field 
gradients,  and 

iii)  accelerating  electric  fields  in  the  center  that  are 
enhanced  by  radial  focussing. 

The  qualitative  description  given  above  is  substantiated  by  the 
analyses  in  the  appendices.  Here  we  give  briefly  their  conclusions. 
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Appendix  I  derives  the  limit  to  the  parallel  electric  field 
9 

E ||  =  10  V/ra  ,  by  balancing  the  Joule  heating  due  to  skin  currents  with 
the  conductive  heat  transfer  into  the  cavity  walls. 


The  exact  expressions  for  the  accelerator  cavity  modes  are  given 
in  Appendix  II.  The  modes  consist  of  a  single  magnetic  field  component, 

(r,z)  ,  which  is  a  product  of  the  axial  and  spheroidal  wave  function. 

The  electric  field  has  an  axial  and  radial  component  given  by  the  curl  of 

(times  c  /iw)  .  Exact  expressions  can  then  be  written  down  for  the 


quantities  of  interest,  e.g. 
total  energy  in  the  cavity,  etc. 


EJ.,wall/Ell,axis 


,  the  cavity  loading,  the 


The  other  appendices  work  out  some  of  the  details  in  the 
formulation,  and/or  contain  numerical  evaluations.  Appendix  III  gives  the 
field  ratio  in  terms  of  axial  mode  number,  and  cavity  frequency 

and  shape. 


Appendix  IV  derives  the  electron  acceleration,  and  does  the 

emittance  calculation.  Here  the  fields  are  approximated  by  their  Gaussian- 

Hermite  form.  The  energy  gain  is  again  a  Gaussian-Hermite  function  of  the 
!< 

frequency  w  .  The  transverse  momentum  transfer  on  axis  vanishes  by 
symmetry,  but  the  maximum  transverse  momentum  gain,  at  about  X/4  off- 
axis,  is  of  the  same  order  as  the  parallel  gain:  therefore  the  radial 
bunch  width  should  be  small  compared  to  a  wavelenth. 

Finally,  Appendix  V  is  a  compendium  of  mathematical  formulae  that 
were  found  useful  in  the  course  of  the  analysis.  One  item  may  be  interesting, 
namely  a  numerical  comparison  of  the  spheroidal  functions  with  the  various 
approximations  used  throughout.'  Appendix  VI  is  a  collection  of  miscellaneous 
physics  results,  and  Appendix  VII  treats  open  cavities. 
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The  overall  conclusion  from  these  analyses,  and  especially  from 
numerical  evaluations  at  MRC  in  Appendix  VIII  (AMRC-N-1S6)  is  that  the 
transverse  RF  accelerator's  performance  is  limited  by  the  ratio 

AE/eE^  L  <  0.5  -  1.0  .  (1) 

Here  AE  is  the  energy  transfer  to  the  electrons,  L  the  cavity  length,  and 
Ej^  a  typical  normal  wall  field.  This  value  is  comparable  to  that  obtained 
in  conventional  cavities,  and  there  appears  to  be  no  advantage  in  using 
cavities  generated  by  extending  optical  resonators  to  the  microwave  regime. 

Although  we  know  of  no  fundamental  limitations  indicating  that 
Equation  (1)  is  a  fundamental  limit  for  all  cavities,  possibly  with  some 
factor  of  order  unity,  it  appears  for  the  moment  that  the  decisive  way  to 
improve  the  accelerating  gradient  is  to  increase  the  breakdown  limit  to  E^. 

No  detailed  analysis  has  yet  been  made  for  purely  propagating 
modes  that  are  geometrically  focussed  down  to  a  small  focal  region.  This 
way  one  has  a  handle  on  the  spatial  field  pattern  on  the  wall,  through 
the  placement  of  the  emitters  on  the  wall.  Therefore,  the  field  direction, 
Ej^/E||,  (with  E ||  at  A/4  away  from  the  wall)  is  under  control. 

Using  a  forced  mode  instead  of  a  cavity  eigenmode  allows  one  to 
select  a  frequency  relatively  independent  from  the  field  pattern  on  the 
wall.  Therefore,  the  coupling  to  the  particles  can  be  somewhat  influenced 
because  the  coupling  partially  determined  by  the  field  frequency. 

If  the  search  for  a  large  coupling  together  with  small  normal 
wall  fields  continues  to  give  negative  results,  we  must  conclude  that  the 
desired  increase  in  accelerating  gradient  should  be  sought  elsewhere,  for 
example  by  increasing  the  breakdown  limit  by  more  suitable  wall  materials, 
or  local  magnetic  insulation. - 
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ORIGIN  OF  RESONATOR  CAVITY 


b.  Cylindrically  Symmetric  Resonator 


c.  Accelerator  Cavity 


FIGURE  I 
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From  A.  E.  Siegman 
"An  Introduction  to 
Lasers  and  Masers" 


HG.  6*23  Example*  of  higher  ord^ 
Hermite  gaustian  optical  resonator 
mode*  at  photographed  m  the  output 
beam  of  a  later  oscillator. 

(linear) 


Figure  2 
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Accelerating 


a. 


GJt  =  7T+  € 

Accelerating 


c. 


Figure  3. 
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APPENDIX  I 

TRANSVERSE  RF  PEAK  FIELD  STRENGTH  ESTIMATES 


The  peak  electric  field  strength  E  normal  to  a  conducting  surface 

7 

ls  typically  on  the  order  of  10  volts  per  meter.  This  value  is  set  by 
breakdown  due  to  electron  emission  and/or  particle  bombardment.  Larger  E^'s 
are  possible,  for  example,  by  magnetically  insulating  the  conducting  surface 
in  the  high  gradient  region  with  a  parallel  magnetic  field  perpendicular  to 
the  electric  field.  The  accelerating  gradient  of  particle  accelerators  using 
normal  fields  is  comparable  to  eEjX:  any  technology  that  increases  the 
attainable  E^'s  could  have  an  impact  on  the  maximum  acceleration. 

However,  it  is  also  possible  to  produce  an  accelerator  with  high 
magnetic  fields  strengths  B|  parallel  to  the  wall,  and  negligible 
normal  electric  field  components.  This  may  be  accomplished  by  utilizing 
mirrors  to  reflect  a  travelling  wave,  as  is  common  practice  in  laser  work. 

The  wave  direction  is  normal  to  the  wall,  and  the  magnetic  field  (and  the 
accompanying  electric  field)  is  perpendicular  to  the  propagation  direction. 

Normal  reflection  from  a  good  conductor  results  in  a  very  small 
electric  field  at  the  surface,  since  the  reflected  wave  interferes  destruc¬ 
tively  with  the  incoming  wave  at  the  conducting  surface.  The  failure 
mechanism  for  such  a  surface  is  different  from  that  of  electric 
breakdown.  For  the  regime  of  interest,  failure  is  due  to  the  surface  layer 
heating  by  the  currents  induced  by  the  incoming  wave.  This  is  known 
as  the  skin  depth  effect,  since  the  electric  and  magnetic  fields  only  pene¬ 
trate  approximately  one  skin  depth.  This  failure  mechanism  scales  with 
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frequency  and  conductivity  of  the  reflector;  it  also  depends  on  the  heat 
transfer  rate  from  the  skin  depth  layer  to  the  interior. 


The  reflecting  power  for  normal  incidence  on  a  plane  conducting 
surface  is  R  =  1  -  S  ;  Here  S  is,  approximately. 


/l  6ti 

s  '  V— ; 


e  v 
o 


where  v  is  the  frequency  (in  Hz)  and  a  the  electrical  conductivity 

(in  (ftm)  *).  The  energy  per  surface  area  absorbed  in  the  conductor  is  equal 

1  2 

to  S,  times  the  average  electric  density  ^€Q|Ep|  (Ep  is  the  peak 
field  strength,  in  V/m) ,  and  the  pulse  length  T  times  the  light 
speed  c. 


The  conductivity  of  copper  is  function  of  temperature  is  shown  in 
Figure  1.  The  conductivity  changes  one  order  of  magnitude  from  room  tempera¬ 
ture  T  =  300°K  to  the  melting  point  at  about  1000°K.  A  reasonable  average 
value  over  this  temperature  range  is  o  -  3  x  10^  (ftm)  *.  The  conductivity 
could  be  increased  by  an  order  of  magnitude  when  the  mirrors  are  cold,  10°K 
say,  at  the  beginning  of  the  electric  field  buildup. 


The  energy  is  deposited  as  heat  over  one  electrical  skin  depth. 


the  heat  diffuses  inward  over  a  characteristic  distance 


At  /ire 
P 


v 


where  A  is  the  thermal  conductivity,  and  cy  is  the  heat  capacity  at 
constant  volume.  For  copper  at  room  temperature  the  thermal  conductivity  is 
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A  =  40 0  (Jn"1sec’1*K~^) ,  the  heat  capacity  cy  =  3.S  x  106  (Jm  ^  °K  *). 

For  oscillating  heat  load  a  "thermal  skin  depth"  can  be  defined  in  analogy 
to  the  electrical  skin  depth.  Setting  tn  -  2ir/tp  , 

V5th =  : 

the  numerical  value  for  copper  at  room  temperature  is  about  10. 

Since  the  thermal  skin  depth  is  an  order  of  magnitude  smaller  than 
the  electrical  skin  depth,  there  will  be  little  influence  from  heat  diffusion 
over  a  typical  oscillation  timescale  T  =  2tt/oi -  After  about  100  oscilla¬ 
tions,  however,  the  heat  conduction  starts  to  come  in,  and  for  even  longer 

times  the  heat  conduction  dominates.  The  frequency  and  pulse  width  envisioned 
for  the  RF  cavity  is  about  1GHz  and  Tp  -  10  ^  sec,  or  about  1000  oscilla¬ 
tions.  Hence  it  is  marginally  reasonable  to  assume  that  heat  conduction  can 
spread  the  heat  over  the  thermal  diffusion  length  6^  =  /h/vc^  Tp  . 

Mirror  failure  is  defined  by  the  melting  of  the  reflector.  This 
occurs  when  the  input  energy  has  increased  the  mirror  temperature  to  the 
melting  point,  about  105oK  for  copper.  Therefore,  the  peak  field  is 
limited  by 


j-Eo  -Up! 


where  AT  is  the  temperature  increase.  Inserting  the  numerical  values  leads 
to 


Ep  =  3  X  109  (Tpv)"*  (V/m)  . 

Some  limitations  of  the  above  analysis  are  obvious:  the  neglect  of 

the  spatial  dependence  exp  -  2 x/6  in  the  heat  generation,  the  assumption 

0 


that  the  heat  is  spread  out  evenly  over  a  thermal  diffusion  length,  and  the 
decrease  of  conductivity  with  increasing  temperature  all  tend  to  increase  the 
peak  field  estimate.  The  assumption  that  the  peak  field  is  attained  instan¬ 
taneously  instead  of  increasing  linearly  with  time  decreases  the  estimate. 

A  working  estimate  for  our  purposes  will  be 

Ep  =  109  (Tpv)"*  V/m  . 

Data  from  laser  mirror  failure  tend  to  corroborate  the  above  estimate. 


A  more  complete  analysis  would  use  the  heat  diffusion  equation 

'  2  2 

S  If  '  sm  = 


analytic  solution  of  this  equation  is  not  possible  except  under  simplifying 
assumptions3  (a  constant,  etc)  similar  to  the  ones  made  above. 


oxlmate  behavior  of  electrical  conductivity 
copper  as  a  function  of  temperature  (From  Ref.  2) 
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APPENDIX  II 


ELECTROMAGNETIC  FIELDS  IN  ELLIPSOIDAL  CAVITY 


The  electromagnetic  fields  envisioned  in  the  transverse  RF 
accelerator  have  the  electric  vector  mainly  parallel  to  the  cavity  axis, 
while  the  magnetic  field  is  roughly  perpendicular  to  the  axis.  These 
conditions  are  satisfied  by  one  case  that  is  exactly  solvable,  namely 
the  purely  solenoidal  oscillations  mentioned  by  Meixner,  (Reference  1, 
Sections  4.41  and  4.43).  The  particular  case  where  all  fields  can  be 
derived  from  one  magnetic  field  component  (see  Figure  1)  will  now 

be  discussed. 


Maxwell's  equations  in  vacuum  are 

V  x  E  =  - 
v  x  t  3t  , 

(la) 

(lb) 

here  is  the  current  density,  and  all  quantities  are  in 

MKS  units.  It 

is  simpler  to  work  with  a  system  where  various  factors  c 
defining. 

are  absorbed: 

B'  =  c B.  ,  J  =  poci  , 

(2) 

and  dropping  the  primes,  equation  (1)  becomes 

OxE.-I  f&  .  " 

(3a) 

—  c  3t 

(3b) 
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These  equations  can  be  combined  into 

■  «> 

The  free  oscillations  of  the  ellipsoidal  resonator  are  obtained 
by  setting  J  =  0  and  assuming  a  harmonic  time  dependence  exp  -  iu>t. 

V  x  7  x  B  =  k2B  .  (5) 

with  k  =  Cl)/ c ,  and  with  the  boundary  condition  that  the  electric  field 
parallel  to  the  cavity  wall  be  zero. 


In  an  arbitrary  orthogonal  curvilinear  coordinate  system, 
equation  (5)  for  the  Bs  component  of  the  magnetic  field  becomes: 


h,h2  k2Bj 


♦ 


3  ha  3  h  o 
3qi  h3h,  3q,  n3B* 


+ 


_3  h,  3 
3q2  h2hj  3q2 


h,Bj 


0. 


(6) 


For  prolate  ellipsoidal  coordinates  (see  Figure  2),  where 


qi  -  n.  qa  -  C,  qs  *  ♦  ; 

x  *  r  cos  \l>,  y  =»  r  sin 


r  »  |  (l-n2)**  (CM)5* 

*  *  §  n  £ 


(7) 


where  f  is  the  focal  distance,  the  coordinates  and  scale  factors  are 

\  -  f  Cc*-n2)%  /  <C2-D%  ; 

h*  -  |  (C2-n2)H  /  (l-n2^  , 

h,  »  |  (l-n^tt2-!)5* 


3-2 


Then  equation  (6)  becomes 


(C2-n2)  Bj*  (l-n2)1*  (l-n2)^,*  (C2-!)*5  ^  (C2-1)5sb3  =  o. 


C9) 


which  is  clearly  separable,  but  a  non-standard  form  of  the  spheroidal 
equation.  Kith  the  relation  (equation  for  £  ) 


Bs 


[(CM)15  ^  -  [jf  (EM)  j|  - 

one  obtains  the  spheroidal  wave  equation  for  m=l  : 

E3  trl  3  3  2,  5  fc2  2.  /kf \2  1 

3£  ^  *  3n  ^1_n  ^  3n  +  ^  'n  M  2/  ’  C2^! 


(10) 


B3  =  0.  (11) 


Hence  the  exact  solution  is  given  as  a  product  of  two  spheroidal  wave 
functions  and  their  derivatives: 


B^C^jil)  «  const.  (h,£)  Sj^  (h,t|)  ,  (12) 

where  h  is  the  parameter 

h  =  ^  .  (13) 

The  electric  field  component  parallel  to  the  cavity  wall 
£  *  3  constant,  is  : 

:  (14) 

the  boundary  Condition  is  that  Ej)  *  0,  hence  a  somewhat  complicated 
condition  on  the  zeroes  of  the  first  derivative  of  the  radial  spheriodal 
function,  as  will  be  seen  below.  Numerical  work  appears  to  be  necessary 
for  firm  results.  The  normal  electric  field  is 
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(15) 


h  s  h  v/  £2-n2 


3 

3n 


rl-n  B. 


These  formulas  give  the  exact  fields  in  an  ellipsoidal  cavity. 


An  attractive  feature  in  the  transverse  RF  accelerator 
is  the  reduction  of  normal  electric  field  at  the  wall  £  =  compared 
to  the  maximum  field  strength  that  accelerates  the  particles  along  the 
z-axis  at  5*1*  Therefore  one  should  evaluate  the  maximum 


Ee|s-g0  - 

'n  E=1  max.n 


Rln(-h’?o) 

(1 

-n2)^  Sln(h,n)j 

max 

2 

(l-n  )  sln(h.n) 

3 

3E 

max 

>/e2- 1  Rln(h,S) 

^o 

,(16) 


in  addition  to  the  energy  transfer  to  accelerated  particles  given, 
essentially,  by  the  Fourier  transform  of  equation  14. 


Before  continuing  with  the  exact  formulas  (14)-(16),  a  few 
remarks  are  in  order. 


i.  The  derivation  of  equation  11  is  not  elegant:  another 
derivation  is  as  follows..  Consider  the  relation  (problem  13.7  of 
Reference  6) 

-7  x  V  x  =  e^  jV*  -  •  (17) 

where  e^  is  the  unit  vector  in  the  azimuthal  direction,  and  'j,®4(r,z)  is 
rotational ly  symmetric.  (This  relation  is  easily  verified  in  cylindrical 
coordinates)..  The  -<f»/r2  term  can  be  conceptually  generated  from  the 
azimuthal  mode  cos  $  (m=l)  in  the  Laplacian.  Hence  equation  (11) 

follows  directly. 
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ii.  The  functions  U  (n)  *  (l-n2!5  S.  (h,n)  and 

.  n  in 

Vn(5)  =  (£*-1)^  Rlnfh,C)  are  entire  functions:  in  particular  near 

£,n  =  1  there  are  no  singularities.  For  numerical  purposes  it  may  be 
nicer  to  integrate  the  differential  equations  for  Un  and  Vn  in  a 
search  for  the  correct  boundard  condition  E^  =  0.  The  appropriate 
equation  is  a  separated  version  of  equation  9. 

a2u 

(1-H2)  - n  -  (A2-h2n2)U  =  0  ,  (18a) 

3n2  n 

and  the  same  equation  for  Vn(£).  The  boundary  condition  E^  =  0  is 
now 

=  0  .  (18b) 

«o 

The  eigenvalue  Ar  =  A^  is  the  one  from  the  usual  spheroidal 
equation,  and  they  are  known,  (eg.  Reference  3,  Page  29,  reproduced  in 
Table  1).  The  resonance  frequencies  are  probably  given  in  Reference  4, 
but  this  reference  is  not  easily  assessible. 

The  simpler  estimates  performed  in  the  initial  stages  of  the 
field  and  energy  transfer  calculations,  (eg.  Reference  5)  compare 
favorably  to  the  exact  electromagnetic  fields  given  above,  and  to  the 
energy  transfer  calculations  that  follow.  First,  let  us  determine  the 
frequency  from  the  boundary  condition,  *  0  ;  with  equation  (14), 

ri„  0’»]  -  o  ;  <i9) 

in  a  future  numerical  calculation  we  will  solve  for  the  frequency  parameter 
h  »  (£0)  as  function  of  axial  mode  number  n  and  cavity  wall 

location  CD*  Here  the  WKB  approximation  given  in  Reference  1,  Section  3.91, 

will  be  used.  Approximately, 


Rln  <»•«  *  {if  {if  J.  to  •  (J°0 

where  Jj  is  the  Bessel  function  of  first  order, 

s  =  s(h,n;5)  =  \  arc  cos  ^  "  ^f7- '  (20b) 

lc 

and  q  =  2n-l.  The  factor  (s/s1  )”*  is  a  slowly  varying  function  of 
5.  viz. 


1  - 


.qjLV 

16h^C 


1  - 


_a__ 

ZhE? 


-A 

2h 


arc  cos  f 


-1 


\ZV^l 

(q**3)  (  e~2 T 

16h*£J 


(21) 


for  £»1  this  becomes  approximately  s/s'  =  £.  Therefore  it  is 
reasonable  to  ignore  this  factor  in  (20a)  .  The  boundary  condition  En  =  0 
then  yields  the  result 


s 


(22) 


or  approximately,  neglecting  the  last  term 

hi  “ h  •  (*0-1)*  (ju+  arc  cos  e;1)  '•  (23a) 


Not  neglecting  the  last  term  gives  h^  »  !s(h1^\hj-*-(q2*3)/4^2) 

Here  l  »  1  corresponds  to  a' half-wave  cavity,  and  1=2  to  a  3/2-wave 
cavity. 
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The  exact  electric  field  along  the  axis  is 


En  =  Ep  Nn  (l-n2)_!i  Sln(h,n)  ,  (24) 

where  E  is  the  peak  field  value,  and  N  is  a  normalization  constant. 
P  n 

-U 

Asymptotically,  for  h»l,  in  a  region  of  width  h'^  about  zero7 


E  =  E  N  H  (n/fr)  exp  -  h  r\2/2  ,  (25) 

n  p  n  r  v 

where  r  =  n-1  (in  general,  S  —  H  ) 
v  s  ’  mn  n-nr 

The  previous  estimates  of  the  energy  transfer  were  based  on 
equation  (25),  but  with  r  =  n;  however,  in  that  case  the  frequency 
estimate,  hg,  was  given  by 

he  *  [jo*  +  ^T-  arC  COS  C]  *  (26a) 

Hence  for  a  corresponding  field  pattern,  n  in  equation  (26)  should  be  re¬ 
placed  by  n-1;  then  the  difference  between  equations  (23)  and  (26a)  is 

4hi  *  id  -  :  t26b) 

the  first  two  values  are  Ahj  =  0.56  and  Ah2  =  0.19,  for  a  relative 

change  in  frequency  Ah/h  less  than  about  10%.  The  conclusion  here  is 
that  the  energy  transfer  as  computed  previously  is  consistent  with  a  parallel 
compulation  using  the  more  accurate  vector  electric  fields. 

A  conclusive  energy  transfer  calculation  must  be  done  with  accurate 
numerical  data  on  the  spheroidal  functions,  and  their  zeroes.  These 
calculations  are  underway,  but  for  the  moment  it  is  reasonable  to  trust  the 
previous  energy  transfer  rates.  It  should  be  mentioned,  however,  that  the 
approximation  to  the  axial  spheroidal  function  S  j  ^  Hn_^  exp  ...  is  not 
strictly  valid7  in  the  region  of  exponential  decrease  [ >>1  *  and  this 
fact  could  effect  the  energy  transfer  calculation. 


Ej^ at  the  cavity  wall. 


An  important  quantity  is  the  normal  field  at  the  cavity  wall. 

Elsewhere  it  is  shown  that  the  parallel  electric  field,  or  rather  the 

g 

equivalent  magnetic  field  at  the  wall,  can  be  up  to  10  V/m,  while 
electrical  breakdown  sets  in  at  30  times  lower  values,  approximately 
3  *  107  V/m.  Therefore  it  is  important  that  E^/E^/  be  small. 

The  ratio  between  the  normal  field  at  the  wall  and  the  parallel 
field  along  the  axis  is  given  by  equation  (16).  This  formula  is  being 
evaluated,  and  the  results  will  be  given  in  another  appendix:  here  an 
estimate  of  Ej //  will  be  given  with  simple  physical  considerations. 
Namely  the  existence  of  a  spatial  magnetic  field  pattern  on  the  cavity 
wall  leads  to  a  comparable  surface  current  pattern,  hence  to  a  surface 
charge  density  which  finally  generates  the  normal  field. 

Figure  3  gives  the  accelerating  cavity  with  the  directions  of 
the  main  field  components  E yy  along  the  axis,  and  in  the 

midplane  parallel  to  the  cavity  wall.  The  surface  current  J  induced 
in  the  wall  is  approximately  given  by 

%  =  V  *  (27) 

found  by  integrating  the  magnetic  field  along  the  sides  of  the  rectangle 
given  in  Figure  4a.  The  spatial  dependence  of  the  =  B^(n)  gives  the 

charge  accumulation  Q  in  the  cube  of  Figure  4b: 

J  =  1//  £*lii  •  t28) 

here  l ^  is  a  typical  length  scale  in  the  variation  of  B^, 

O  is  the  surface  charge,  and  w"1  is  a  typical  time  scale.  The 
surface  charge  gives  the  normal  electric  field  according  to 


El*^ 

o 
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The  normal  electric  field  at  the  cavity  wall  EA  ,  compared  to  the 
peak  parallel  electric  field  \/U  away  from  the  wail,  Ej  =  c  , 

Ea/e,  =  («-sk)  '* 

where  k  =  m/c  is  the  free- space  wavenumber. 

The  above  estimate  shows  that  the  large  fields  predicted  for  the 
transverse  RF  accelerator  are  obtainable  when  the  free-space  wavelength 
is  small  compared  to  the  wavelength  of  the  fields  on  the  wall. 

The  accelerating  field  along  the  axis  in  the  cavity  center  is 

the  parallel  field  at  the  wall  times  an  additional  radial  factor  4R/A  , 

where  R  is  the  cavity  radius  and  A  is  the  wavelength.  The 

magnetic  field  B .  goes  with  radius  r  roughly  as  1/r  because 
f  ♦ 

7— ■—  =  271  r  (r3  is  a  constant,  for  each  integration  path  along  the 

maximum  of  B.  .  The  electric  field  in  the  center  is  roughly  equal  to 
0 

cB.  along  the  smallest  integration  path  encircling  the  center,  at 
0 

r  =  A/4  .  Therefore  the  normal  field  at  the  wall  E^  is,  compared 

to  the  accelerating  field  E,  on  axis, 

Ej/EH  =  /8 jj R  .  (31) 

Therefore  the  breakdown  limitation  can  be  avoided  by  using  wavelengths 
short  compared  to  the  cavity  dimensions.  An  estimate  based  on  the  exact 
formula  (16),  gives  elesewhere,  is  in  general  agreement  with  (31). 
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CAVITY  LOADING 


The  electron  beam  that  is  accelerated  by  a  particular  externally 
excited  cavity  mode  can  itself  excite  other  cavity  modes.  These  modes,  once 
built  up  to  large  amplitudes,  have  their  own,  probably  deleterious, 
influence  on  the  accelerating  beam.  The  mode  excitation  due  to  a  prescribed 
electron  beam  is  given  in  the  sequel.  _ 

Various  important  simplications  are  made.  The  beam  electrons 
are  relativistic,  and  therefore  they  keep  their  velocity  approximately 
constant,  equal  to  the  light  speed  c.  This  approximation  necessarily  ignores 
any  radiation  from  beam  acceleration;  most  of  the  radiation  is  peaked  in  the 
forward  direction,  and  will  escape  from  the  cavity.  It  will  be  treated  in  a 
separate  calculation. 


The  beam  is  supposed  to  be  injected  into  the  cavity  along  the 
z-axis,  preserving  the  azimuthal  symmetry  already  imposed  on  the 
accelerating  mode.  Then  the  magnetic  field  generated  by  the  beam  current 
will  be  in  the  same  direction  as  the  cavity  mode  field,  namely  the 
azimuthal  direction. 

The  above  considerations  allow  simplification  of  the  full 
inhomogeneous  vector  wave  equation  to  an  inhomogeneous  scalar  wave  equation 
for  the  magnetic  field  *  e^  b  (See  equation  17  sqq.) 

72b  -  P-  b  -  0  3  C.  (32) 
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is  the  unit  vector  in  the  azimuthal  direction,  and  C  is 
component  of  the  curl  of  the  current  source  term, 

e^  COSt)  =  Po  c  7  k  ^  j(*,t)  (33) 

where  j(>St)  is  the  current  density.  A  reasonable  functional  form 

for  the  n'th  electron  bunch  is 

jn(X,t)  =  exp  ['  (Vct)2/2Lz  ‘  ]  <34) 

where  I  is  the  peak  current,  Lr  is  a  typical  width  in  the  radial 

direction,  z  is  the  center  of  the  particle  bunch  along  the  cavity 
n 

axis,  with  a  typical  bunch  length  L  .  It  is  easy  to  verify  that 
Vxj_  's  only  component  indeed  points  in  the  e^  -  direction. 

A  suitable  set  ^  of  real  eigenfunctions  for  the  cavity  is 
n 

found  above.  The  inhomogeneous  equation  (32)  is  then  solved  by 

b  (*,t)  =  /  Gfx.t;^,^)  CCx,.^)  ,  (34) 

where  the  Green's  function  G  can  be  expanded  in  terms  of  the  eigen¬ 
functions  (e.g.  Reference  6,  equation  7.3.24) 

b  ■  c2  ?  Yj(x)/“j  /  sin  jV.(t-tQ)  ]  JdVj(V  C(W  *  (35) 

Here  the  summation  is  over  all  indices  of  the  eigenfunctions,  for  the  cavity 

j  *  (n.A) 
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is  excited 


These  expressions  basically  say  that  each  mode 
by  the  source  C  with  a  coupling  strength  that  depends  on  the  spatial 
overlap  of  the  mode  with  the  source,  and  on  the  temporal  overlap.  All 
that  is  left  to  do  is  to  evaluate  these  overlap  integrals  explicitly. 

When  this  is  done  one  can  go  back  and  consider  the  effect  of  these  beam- 
excited  modes  on  the  beam  itself. 

In  the  evaluation  are  needed  the  explicit  forms  of  the 
eigenfunctions,  the  current,  and  the  frequencies.  The  eigenfunctions  are 
given  in  equation  (11): 

Vi>  ’  ""I  Sln  Rln  <h„£  «  •  <36> 

where  A  .  is  the  normalization  constant. 


h  is  the  eigenvalue,  related  to  the  frequency  by  h  * 

indices  (h  =  h.  „)  are  omitted.  The  eigenvalue  satisfies 
lnx 


its 

(37) 


The  source  C(xQ,to)  i-n  spheroidal  coordinates  is 


(38) 


Notice  that  this  formula  is  for  one  electron  bunch  only:  for 
there  is  another  summation  over  bunches,  indexed  by  j,  say: 


N  bunches 

C  *2.  C.  . 

J  J 
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With  these  expressions  it  is,  in  principle,  possible  to  evaluate 
equation  (35);  howeyer,  the  calculations  are  terribly  messy  and  they  are 
not  performed  now.  Instead,  we  make  the  following  physical  arguments  in 
an  attempt  to  search  an  acceptable  conclusion. 

The  first  observation  is  that  the  spatial  overlap  integral  can 
be  approximated  in  two  limits,  one  in  which  a  typical  length  scale  of  the 
eigenfunctions  is  much  larger  than  the  scale  of  the  charge  bunch;  the 
other  limit  applies  when  the  opposite  is  true.  These  are  the  low  and  high 
mode  number  cases,  respectively.  In  the  first  case  the  integral  is 
approximately  constant,  in  the  second  case- the  integral  is  zero.  The 
upshot  here  is  that  only  a  finite  number  of  modes  is  coupled  to  the 
source  and  that  the  coupling  constants  are  all  of  the  same  order.  In 
this  connection,  notice  that  the  normalization  constant  is  of  order 
(peak  value  of  the  function  x  integration  interval) 2/j,  hence  that 
the  normalization  and  peak  value  cancel  each  other,  except  perhaps  for 
some  dependence  on  the  integration  interval  in  the  normalization. 


It  is  a  different  story  for  the  time  integration.  Here  the 
coupling  involves,  mainly, 

/dt<> sin  [»«*  <*-*<>>]  “p  -  (i  vct)’/,Lz 


(39) 


which  can  be  large  only  when  the  correct  bunch  spacings 

relative  to  the  frequency  of  the  accelerating  mode  are  obtained.  Otherwise, 
the  integrand  in  (39)  will  be  oscillatory,  and  vanish  when  integrated  over 
many  bunches:  one  bunch  may  excite  a  mode  to  a  certain,  small,  amplitude, 
but  the  next  bunch,  or  a  bunch  a  few  further  down  the  line,  will  take  the 
energy  from  this  mode:  the  net  result  is  little  excitation  of  extraneous 
modes. 
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The  above  scenario  is  reasonable  if  the  frequencies  are 
incommensurate.  This  is  highly  likely,  but  untested  at  the  moment.  The 
best  estimate  gives  for  the  frequency  a  linear  expression  in  terms  of 
Bessel  functions  zeroes:  these  are  certainly  incommensurate. 

Our  conclusion  is  that  modes  other  than  the  accelerating  mode 
can  be  excited  in  the  cavity,  but  that  these  modes  remain  at  a  relatively 
low  level.  Therefore  they  will  have  a  minor  effect  on  particle  acceleration 
in  the  transverse  RF  accelerator. 
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Further  Analytical  Estimates. 

The  total  electromagnetic  energy  W  in  the  cavity  can  be 
easily  estimated  as 

W  =  E2V 
2  P 

times  a  geometrical  factor;  Ep  is  the  peak  field,  V  is  the  volume. 

The  geometrical  factor  of  order  unity  must  be  evaluated  using  the  particulars 
of  each  cavity.  Remember  that  in  radiation  field  the  time-averaged 
electric  and  magnetic  field  energies  are-equal, 

«0  / |E|a  dv  =  p"1/!!)2  dv  .  (40) 

and  therefore  it  is  sufficient  to  evaluate  the  least  cumbersome  of  these 
expressions. 

The  maximum  magnetic  field  energy,  equal  to  the  total 
E-M  energy  W,  becomes  in  spheroidal  coorinates: 

S  r 

W  -  (f)  W0  *1  Np  / dn  Sln  (h’n)  Rl2n  (41) 

1  -1 

Here  Bp  is  the  peak  magnetic  field,  and  the  double  integral  with 
henceforth  be  denoted  by  I.  Np  is  a  normalization  constant  to  be 
discussed  later  Equation  (41)  clearly  displays  the  expected  form,  an 
energy  density  Bp 2/a/z0  times  a  volume  (for  a  confocal  cavity,  the 

diameter  is  equal  to  the  focal  distance  f,  hence  in  this  case 
V  *  irf*),  times  a  geometrical  factor  given  by  the  remaining 
factors . 

The  double  integral  can  be  split  in  two  terms  with 
C*-n*  •  (C*-U  ♦  (i-na)  : 
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I  •  J45  (t!-i)  Rjn  (  hntj5  )  fir,  (i-n!)  s’ln  (h.n) 

1  -1 

4  1 

JdCRln  C  hn£jC  }  /dn  C1'n2)  Sln  (h,n)  ’  (42) 


1  -1 

It  is  convenient  to  choose  the  normalization 

1 


/dn  s2n  (h.n)  =  l 

-1 


(43) 


Note  that  this  normalization  differs  from  the  one  used  by  others,  eg. 

f dn  .  .  *  - as  used  by  Meixner  and  Schafke  : 

J  (n+4) (n-l)n 

these  notational  difficulties  are  a  recurring  nuisance.  The  remaining 
n  -  integral  is 

r 

_/dn  d-n1)  s’n(h,n)  »  5^  .  («) 

where  An  is  the  eigenvalue  in  the  spheroidal  differential  equation 


dir  (1-nZ)  an +  q  (A*h2»n;m) 


s  (h.n)  *  0 

mn 


(45a) 


where 

q  »  -  An  ♦  h*  (1-H2)  -  .  (45b) 

Series  for  Ar  in  powers  of  h2  are  well-known  and 
this  fact  makes  relation  (44)  useful.  Its  proof  is  simple:  multiply 
equation  (45a)  for  Ijq  with  (hj.n)  belonging  to  h2  =  h2+Ah2, 

interchange  the  subscript  0  and  1,  subtract  and  integrate  over  the 
interval  1 ,  l"l  .  The  result  is 
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0  =  (Ah2)  fin  jq(A0,h2,n;m)  -  q(Ai  ,h* ,n;m)j  s^Oyn)  »  (46) 

Taking  the  limit  Ah2—- o  gives 

/*>  Ur  SL  <»•'”  -  0  •  (47a) 

which  is  equation  44  when  using  the  explicit  form  for  q,  equation  45b. 

The  eigenvalue  A  in  Table  1  is  related  to  A  by 
6  mn 

A  =  h2-Xran  ;  hence  equation  (44)  becomes 

y\  3A 

dn  (i-n2)  (h»rj)  =  l-3~  •  C*8) 

-1 

In  the  limit  h2— oo  ,  when  S2  (h,n)  goes  to  a  6  -function 

mn 

(or  derivatives  thereof) ,  equation  (48)  becomes  unity,  in  agreement  with 
Table  1. 


The  integrals  over  the  radial  coordinates 
relation  of  the  same  form,  but  with  an  extra  term: 


satisfy  a 


(52-l)  Ri„  (h.C)  = 


31  a2 

^h7  353h* 


XA2-1  R 

V  _ mn 


(49a) 


c=c. 


This  last  term  is  related  to  the  non-zero  boundary  condition  on  R^. 
This  term  can  be  written  another  way.  A  change  in  the  parameter 

h  (frequency)  implies  a  different  coordinate  £q  ■  CQ  (h2)  where 
vanishes  (see  equation  19). 


Then  3/3h2  can  be  replaced  by  3£0/3h2  x  3/3^,  and  the 
second  derivative  becomes,  3a/3£2  :  This  in  turn,  can  be  re-written 
with  equation  (10).  The  result  is 
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mn£ 


32  \/ 52-l  ^nn 


3?  (h2)  X  -h2£2 
__£__  mn  o 

3hz 


SW1 


R  (h,r  .) 
mu'-  n  l' 


(49b) 


The  integral  I  (equation  42)  becomes  simply 


I  =  I 


ln£ 


(50a) 


and  the  total  E-M  energy  in  the  cavity  is 


W  =  (f\  —  B2  N 

V2/  yo  P  1 


2  I2 
p  iln£ 


The  final  ingredient  in  this  formula  is  the  normalization  constant 
Np  (and  the  numerical  evaluation  of  equation  49b) . 


(50b) 


defined  by 


The  normalization  constant  Np  for  the  magnetic  field  is 


"♦  =  Bp  Np  Rln  Sl„ 


(51) 


where  Bp  is  the  peak  value  of  the  magnetic  field,  and  Sjn  and  Rjn 
are  normalized  to  unity  according  to  equation  (43),  and  the  corresponding 
equation  for  R.  Given  this  normalization  one  needs  to  find  the  maximum 
of  Sjn  (h,n)  over  |n|  <1  and  of  Rln  (h,5)  over  1<5<5o 


In  the  axial  direction  a  local  maximum  occurs  at  all  zeroes 
z£  of  Sln/dn;  however,  the  absolute  maximum  takes  place  for  the 
largest  of  these,  z^  *  max^  (z£).  This  can  be  easily  argued  from 
WKB- theory .  Hence , 


sm  0».n> 


max 


Sln  (h*zP 


(52a) 
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where  z*  is  the  largest  zero  of 

D 


3n 


o 


In  the  radial  direction  the  maximum  occurs  closest  to  the  center  £=1- 
This  is  evident  from  the  basic  1/r  dependence  of  cylindrical 
coordinates,  and  can  also  be  argued  on  WKB  grounds.  Hence, 


Rln  (h.n) 


max 


Rln  <»•  *n? 


where  r' ,  is  the  first  zero  of 
nl 

dn 


0 


The  normalization  constant  N  is  then 

P 

N„  •  [“in  Ch.ri,)  sln 

The  hooker  in  all  this  is  that  the  spheroidal  functions  are  not  well 
tabulated,  their  zeroes  are  not  well-known,  etc.  These  numerical  data 
need  to  be  calculated  before  one  can  determine  the  numerical  value  of 
the  various  quantities  such  as  equations  (54)  and  (16). 
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(52b) 


(53a) 


(53b) 


(54) 


r 
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preceding  page  #  20  &  21  are  not  film 


TABLE  1 


The  eigenvalues  for  the  spheroidal  wave 

equation  as  function  of  the  parameters  c  =  h  =  fk/2 

M* 


2  + 

«  + 

,2  +  !«• 

20  +  %'• 

»  +  g<* 

«  +  Ti‘' 

“  +.  s- 

22  ♦  p 


4  . 

,  8 

1,029  c 

+  554,631  ** 

152  ^ 

115.568 

111,375  ^ 

977.315,625  ^ 

7,064  „ 

462.736 

5^93+929  ^ 

35,188,194,041 

1.108 

1,245,699  ” 

301.160  ^ 

96,630,117,129 

5.684 

8,485.125  c 

495.016 

487,682,559,375 

1.167,376 
2,255,916,949 
31.984  „ 

78,021,125 


45^  “ 

c* 


132,492,064 
330,543,719.118.327' 
96.594,016 
534,442,755,721,875  ' 


—  124  -c._ 

58.953.125  **' 

5,420  ,  _ 

9,539,098,569  C 
34,094.936 

15,964.450,734,375  c  +  •  ■  ‘ 
339,769.051,816 

599,396.295,493,692.697  ^  +  ’  ' 
69,222,847,348 
2.052,497.609,859,563,685  ^ 
5.894.701.468 

2,302.259,141,535,890,625  ^  “  ‘ 
3,149,302.579.970,160 
5,406,455.2 19, 1 58,600,21 1 ,866,5 
16,316,385,212.272 

23.965.276.291.519.365.828.125 


For  general  m  and  n  :  q  =  2n  -  2m  ♦  1 ,  Y  =  h 
(Meixner  5  Schaflce,  Section  3.2S): 

2h2  -  A  (h)  *  yf +  ***---- [?*+ 5]- 

mn  ° 

-  £5  (f*  +  26 ?’  +  21)  -  384m* (?'  +  1)]  - 

-  pr  [liTTir 133  ?i+ 1 594  **+  5621  ?) "  W  (37?*+ 167?)  +  T  ?1  “ 

-  ji  -2Sf.3if  (639*  +  49409*+  43  3279*  +  22470)  - 

--^(1159*+ 13109* +  735)  +  ^  (9* +1))- 

-  ^[735*“  (S279T  +  61  5299*+10439619»  +  2241  5999)- 
- (57399*  +  12755°9*  +  29895!  9)  + 

+  (3S5f*  +  15059)  -  9)  + 


± _ 1 _ 

y*  11024-1024 

- 

32-1024  1 
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1  <  c  <  », 

-1  <  n  <  1 . 

0  <  4>  <  2n 


Figure  2.  Prolate  Spheroidal  Coordinates. 
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transverse  RF  accelerator  cavity  with  the  irid]or  e  led  roroa^net.  ic  t  ic'id  components 


APPENDIX  III 


Normal  Fields  With  Exact  Mode  Structure 


The  transverse  RF  accelerator  allows,  in  principle, 

extremely  high  accelerating  gradients  because  the  accelerating  electric 

fields  are  mainly  parallel  to  the  cavity  wall.  The  parallel  field  is 

q 

then  limited  by  Joule  heating  to  about  10  V/m  .  This  value  is  much 
higher  than  the  maximum  normal  field  that  can  be  sustained, 

E1  <  J  x  lo'  V/m  .  This  section  evaluates  the  exact  expression  for 
the  normal  field  at  the  wall,  compared  to  the  accelerating  field  along 
the  cavity  axis. 

The  expression  to  be  calculated  is 


(Equation  16  of  "Electromagnetic  fields  in  Ellipsoidal  Cavity").  It 
consists  of  a  radial  part,  with  independent  variable  £  ,  and  an  axial 

part,  with  variable  n  •  Let  us  consider  these  factors  separately,  using 
the  asymptotic  form  of  R  and  S  for  large  frequency  (h  *  ^  »  1) 
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The  radial  eigenmode  R^n  is  approximately 

Rln(h-5’  ’  Jl(s)  •  (2a) 

where  to  lowest  order  in  h 


s  =  h  /s2-l  *  (2b) 

and  Jj  is  the  Bessel  function  of  first  order. 

The  frequency  is  related  to  the  coordinate  £  of  the  wall  by 
the  implicit  relation 


3 

3C 


Rln(h'0| 


(3) 


Using  the  approximation  (2a)  this  condition  takes  place  approximately  when 
vanishes,  or  equivalently,  at  the  maximum  of  .  In  the  region 

s»  this  maximum  is  simply  the  dominant  term  in  the  asymptotic 
expression  of 

J1»M  'So-11'"  •  (4) 

The  additional  factor  is 


£ 

s' 


(5) 
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The  value  of  the  ^-dependent  denominator  in  Equation  (1) 
can  be  found  in  the  opposite  limit  J^(s)  =  s/2  for  s — wO 

Then 


hence 


a  r 

as  M 


C2-l  Rln(h,C)t 
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(6) 


(7) 


The  axial  factor  can  be  evaluated  from  the  WKB  expression  for 
the  function 


yiV  sln(h,n)  *  Hn-i(nv^  «p  -  hn2/*  ;  (8a) 

in  this  expression  one  can  ignore  the  factor  vl-n2  »  at  least  for 

those  large  values  of  h  for  which  (8)  is  valid.  The  differential 
equation  for  (8)  is 

*  Q(n,x)gn  »  0  ,  (8b) 

where  x  a  H\/h~  ,  and  Q(n,x)  *  2n-l-x2 

The  extremes  of  the  function  S^n  lie  between  the  turning 
points;  in  fact,  the  maximum  is  the  local  maximum  closest  to  the  turning 
point.  In  contrast,  the  maximum  of  the  derivative  occurs  in  the  middle 
between  the  turning  points.  This  can  be  shown  easily  from  the  WKB- 
expansion  about  the  turning  point,  x  *  x£  ♦  Ax  : 

g„  =«  /T  a'  Ai(a'/jAx)  »  (9a) 
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"T 


where  a  =  -  3Q/3x  at  the  turning  point,  and  Ai  is  the  Airy 
function.  The  maximum  occurs  at  the  maximum  of  Ai(t)  =  0.5356  ,  at 

t  *  1.003 


The  derivative  is  calculated  with  the  WfCB-expression  between 
turning  points: 


gn(x)  =  Q’%  cos 


[/V 5"dx] 


(9b) 


with  derivative 


g'  =  Q  sin 


in  [M“dx] 


(9c) 


The  maximum  of  this  is  obviously  where  Q  is  largest,  at  x  *  0 


The  maximum  (9a)  becomes,  using  the  turning  point  x*/  2n-l  , 


8, 


n.max 


Vir” Ai  (max) 


(10a) 


while  the  maximum  of  the  derivative  becomes 


®n,  max 


•  (2n-l) 


(10b) 


The  axial  factor  becomes  finally,  using  3/3r»  -  Vh“  3/3x  and 
setting  n  *  0  in  the  square  roots  that  multiply  S, 


In 


[  r  -jj 


max_ 


d-n*)  a  sln(h.n) 


I 

(  max 


=  Vh"- 


(n-h) 


V i 


VStT  Ai(max)  /g*-l 


•  (10c) 


The  ratio  between  normal  field  at  the  wall,  and  accelerating  field 


along  the  axis,  Equation  (1),  becomes  from  (7)  and  (10): 


-7T  -1-1 

2  Ai(max) 


hOl(C-l)* 


The  coefficient 


^Ai  (max)j  =  0.3 


(11) 


The  conclusion  from  this  formula  is  that  the  normal  fields  remain 
small  compared  to  the  accelerating  fields  by  a  combination  of  high 
frequencies  h  =  ,  high  axial  mode  number  n  ,  and  large  wall 

distance  .  However,  none  of  these  factors  is  clearly  dominant, 

although  increasing  the  frequency  is  most  effective. 


The  formula  (11)  agrees  approximately  to  the  previous  estimate. 
Equation  (31)  in  Appendix  I: 


7T  C2  1 

2  u2"  SljjR 


(12) 


When  Equations  (11)  and  (12)  are  written  in  corresponding  quantities  the  dis- 

.  f U) 

crepancies  become  clear.  For  example.  Equation  (12)  becomes,  using  h  =  , 

R  *  4  Je-l  »  and  the  typical  wavelength  on  the  axis  (An-'*) 

Mmes  a  factor  from  the  projection  cn  the  cavity  wall,  v,j|  =  X£rf 


4l*  |-T 


(13) 


These  formulas  are  in  good  agreement;  the  frequency  dependence 

- 1  -  2 
is  h  in  both,  and  at  least  for  £0  large  the  formulas  go  as 
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1/3  U 

but  the  axial  mode  number  dependence  is  n  '  versus  n'*  .  It  should 

be  noted,  however,  that  both  Equations  (11)  and  (13)  contain  various 
approximations  that  might  explain  part  of  the  discrepancy. 

Numerically,  it  appears  that,  for  CQ  =  \/"i~  , 


0.5 


(14) 


only  a  factor  n*^  off  from  (13)  but  different  by  n*^  from  (11). 
The  reason  is  probably  that  the  influence  of  the  axial  mode  number  on  the 
frequency  is  entirely  neglected  in  Equation  (2b). 


The  cavity  parameters  that  determine  the  normal  fields  at  the 
wall  are  the  wavelength  A  compared  to  the  focal  distance  of  the 
cavity  f  ,  h  1  =  V(Trf)  ,  and  the  cavity  shape  as  parameterized  by 
The  axial  mode  number  n  has  a  minor  influence. 
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APPENDIX  IV 


Acceleration  and  Emittance  Calculations 


The  exact  electromagnetic  cavity  modes  discussed  earlier 
determine  the  acceleration  of  relativistic  particles  shot  into  the  cavity. 
Three  terms  in  the  Lorentz  force  that  need  to  be  evaluated  are: 


dP|, 

dt 


dp 

at  a-  E 


B 


(la) 


(lb) 


(lc) 


the  last  term  coming  from  v  *  B./c  =  -  e_B,  .  The  field  should  be 

-0-9  —5  <P 

taken  at  the  actual  position  of  the  particle,  z  =  z(t)  .  The 

magnetic  and  electric  fields  are  given  as  a  product  of  two  spheroidal  wave 
functions  and  their  derivatives: 

B^CC.n)  =  const.  Rln  (h,(i)  Sln  (h,n)  ,  (2a) 


where  h  is  the  parameter 


h 


kf  fu) 
2  "  2c 


(2b) 


The  electric  field  component  E^  ,  roughly  in  the  acceleration 
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direction  along  the  z-axis,  is 


(2c) 


The  field  component  E^  ,  roughly  perpendicular  to  the  z-axis,  is 


Ec  =  h  yT'-n2  3n  yi -n2  b 


0 


(2d) 


These  formulas  give  the  exact  fields  in  an  ellipsoidal  cavity  once  the 
frequency,  or  h  ,  has  been  determined  from  the  radial  boundary  condition 

W  •  0  • 


The  particle  trajectory  is  assumed  to  be  exactly  parallel  to  and 

to  pass  through  a  narrow  region  about  the  z-axis  Ini  <  1  ,  5-1 

The  magnetic  field  is  always  exactly  perpendicular  to  the  particle 

direction.  However,  the  electric  field  components  E  and  E,  are 

n  s 

exactly  parallel  resp.  perpendicular  to  z_  at  n  =  0  only  :  for 
n  /  0  the  E^  field  has  a  small  component  in  the  £-  direction  , 
and  vice  versa.  These  components  will  be  ignored  to  a  first  approximation; 
henceforth  Eu  =  E  ,  and  E,  =  Er 

#  n  15 


The  energy  transfer  in  the  parallel  direction  gives 

the  accelerating  energy  AE  =  cApg  »  while  the  momentum  transfer  in 
in  the  perpendicular  direction,  Ap^  ,  produces  the  emittance.  The 
momentum  transfers  are  the  integrals  over  time  of  Equation  1,  evaluated 
along  the  particle  orbit  z  *  *  ct,  £  *  1  : 


Ap, 


ie  E 

_ EL 

o> 


sm(M) 


(3a) 
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(3b) 


ie  E 


APi,E  3 


co 


4Pi,B 


»  E„  yv 


0) 


Wl(g)  f,-  e-ihn  3  Sln(h> 

)  V  l-n2  3n  /Tv 

/ 


n) 


dn  e"lhn  sln(h,n) 


(3c) 


Here  is  the  peak  electric  field,  and  Ny  >  and 

N-  are  the  normalization  constants  that  connect  the  three  fields: 
They  are 


Hl  <«  ■  \  Vlh’5) 

y 

(4a) 

N±  (5)  -  Nb(£)  =  Nn  Rln(h,C)  ; 

(4b) 

These  formulas  all  have  the  same  normalizing  factor 

|-1 

sin(h.n) 


Nn 


yTn* 


max,  n 


5) 


max,  ^ 


(4c) 


where  Umax.x 
with  respect  to  x 
electric  field  E 

P 


denotes  the  maximum  of  the  function  in  brackets 
.  These  factors  are  chosen  such  that  the  peak 
is  the  single  strength  parameter. 


The  Fourier  transforms  in  Equations  (3)  show  various  interesting 

features:  there  is  a  180  degree  phase  difference  between  Apy  and 

Ap^  ,  due  to  the  factor  -i  resp.  i  ,  and  a  90  degree 

phase  difference  with  Apj^  .  In  addition  there  is  an  extra  phase 

factor  from  the  initial  position-  zQ  *  fi^/2  .  Hence  the  complex 

exponential  exp  -  ihn  becomes  exp  -  ih(n+no)  .  The  phase 

n  should  be  such  that  the  accelerator  maximizes  Ap„  .  Since 
o  « 
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the  functions  Sin(*'»r|)  have  a  well-defined  parity  (odd  for  n  even, 
and  vice  versa)  this  maximization  demands  that 

-i  exp  -  ih(n+h0)  *  cos  hn 

when  S,  is  even  in  n  ,  and  sin  hn  for  S,  odd;  what  is 
in  In 

more,  this  replacement  should  be  done  in  each  of  Equations  (3).  Let  n 
be  odd,  for  the  sake  of  argument.  Then  the  exponential  in  (3a)  and  (3b) 
becomes  cos  hn  ,  but  in  (3c)  sin  hn  -  Now  the  derivative  3/3n 
in  (3b)  of  an  even  function  changes  the  parity  to  odd,  hence,  the  integral 
will  vanish.  Likewise,  the  integral  in  (3c)  is  over  an  even  function 
times  an  odd  factor  sin  hn  ,  and  this  integral  vanishes  also. 

Not  all  electrons  enter  the  fields  with  exactly  the  correct 
phase:  the  momentum  transfer  for  an  electron  with  phase  angle  n^ 
will  be  given  by 


4Pi  <V  ■  I  Vi  I  cos  hnb 

E(V  ■  IVi.E1  sin 

4p1BCnb)  -  -  |AmPiBl  sin 


» 

(5a) 

hnb 

9 

(5b) 

hnb 

1 

(5c) 

an  angle  factor  times  the  absolute  value  of  the  momentum  transfer,  i.e. 
the  integrals  (3)  evaluated  with  the  maximizing  phase. 

Notice  that  these  considerations  pertain  only  to  the  field 
structure  along  n  ,  and  do  not  depend  on  the  particle  distance  from 
the  z-axis:  Its  effect  will  come  in -only  through  the  factor  N(£)  in 
Equation  (3). 
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The  integrals  in  Equation  (3)  are,  to  our  knowledge,  not 
expressible  in  closed  form  except  in  the  limit  h-*-00  .  Then  the 

spheroidal  functions  becomes 

sin(h.n)  ^  Hn_L  (n/h)  exp  -  ho2/2 

and  the  factor  (1-n  )  becomes  simply  unity. 

Whence , 

eE  n„  (C)  r^r 

lAP||l  =  2  J-  -  h~  Hn-l(l/h  3  exp  ‘  h/2  ’  (6a) 

eE  n,  (o  rz: 

•^i.E1  =  -  h\/K“  Hn-1  t/K")  exp  “  h/2  *  W 

and 

eE  N.(Q  r^~ 

*  ~'£  u  ~  Y  h~  Hn-1  exp  "  h/Z  ;  (6c) 

The  magnitude  of  the  transverse  momentum  change  is  basically  proportional 
to  the  accelerator's  energy  gain.  A  small  emittance,  therefore,  depends 
crucially  on  the  smallness  of  the  inital  phase  ,  and  on  the 

transverse  bunch  size  £.  through  the  normalization  factors  of 

D 

Equation  (4) . 

Close  to  the  cavity  axis,  i.e.  £  =  1  ,  Equation  (4)  can  be 

evaluated  with  the  limiting  form  for  the  radial  eigenfunctions 

Rln(h,£)  near  £  *  1  .  These  are  approximately,  for  large  h  , 

»ln<h'E>  -  If?  f  •  (7a) 

and 
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where  the  radius 


In 


(7b) 


The  maximum  in  the  radial  dependence  of  parallel  electric  field 
E^  takes  place  at  the  cavity  axis,  hence 


a 

as 


In 


max,  £ 


(7c) 


The  normalization  factors  (4)  are  near  the  axis: 


sm  Ch.n) 


y  l-n2 


max,  n 


and 


N1  =  Nll  f 


(8a) 


(8b) 


This  gives  for  Equation  6: 


AP 


II, E 


(X/TT)  APB 


(8c) 


It  is  satisfying  that  these  results  are  so  obvious;  in  fact,  the 
equivalent  of  Equation  (8)  has  already  been  used  to  estimate  the  emittance. 
The  present  analysis,  therefore,  adds  the  connection  with  the  spheroidal 
cavity  modes,  but  no  new  insights.  Perpendicular  momentum  transfer  is 
produced  by  the  magnetic  field  ,  but  mainly  by  the  perpendicular 

electron  field  E^  ,  on  account  of  the  additional,  large,  h  in 
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Equation  (6b) .  Taking  only  Ap^  ^  ,  the  perpendicular  momentum  transfer 

is  for  the  edge  of  a  small  bunch  (hn^*  1  ,  hrb<<  1) 

APj_  =  Api)  2ir*  (rb/X)  (zb/\)  .  (9) 

The  conclusion  from  Equation  (9)  is  that  the  electron  bunches 
to  be  accelerated  in  the  RF  cavity  should  be  very  small  compared  to  the 
wavelength  both  in  radius  rb  and  axially  zb  .  Typical 
dimensions  would  be  1%  of  a  wavelength,  in  both  r  and  z  ,  for 

APj/AP||  =  SM  * 
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Energy  Transfer  From  Cavity  Modes 
To  Gaussian-Shaped  Relativistic  Charge  Bunch 


The  energy  transfer  from  cavity  modes  to  a  single  particle  was 
calculated  earlier.  This  section  computes  the  energy  transfer  to  a  Gaussian 
bunch  of  electrons.  The  bunch  is  assumed  to  be  narrow  in  the  radial  direc¬ 
tion,  but  it  has  a  varying  width  along  the  cavity  axis. 

Two  limiting  cases  are  obvious:  a  bunch  that  is  narrow  compared 
to  the  axial  wavelength  of  the  cavity  modes  looks  like  a  single  particle  for 
energy  transfer  purposes.  On  the  other  hand,  a  bunch  that  stretches  over 
a  few  axial  wavelengths  gets  little  energy  transfer. 

Relativistic  motion  of  the  particles  precludes  any  shape  changes 
in  the  bunch.  Therefore,  the  energy  transfer  can  be  calculated  by  averaging 
over  the  particle  density  in  the  bunch, 

n(5)  =  — ^  exp  -  C2/2aV  :  (10) 

aW/2? 

here  N  is  the  number  of  bunch  particles,  C  is  the  bunch  coordinate 
£  ■  z  -  ct,  and  a  is  the  bunch  width  in  terms  of  the  waist  size 


The  energy  transfer  from  cavity  mode  n  to  a  relativistic 
particle  at  a  phase  angle  sQ  *  coAt  *  co^/c  with  respect  to  the  position  of 
maximum  energy  transfer,  is 
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(12a) 


AE(sq) 


AE  cos  s  , 
o  o’ 


where 


AEo 


N 

(eE  X  )  — 
P  n  /2 ft 


l 

/TTW) 

(x  / 

exp  -  j 

(x  / 

(12b) 


The  energy  transfer  to  the  Gaussian  bunch  AE,  is  then 

D 


AE.  =  AE 
b  o 


N 


otW/2? 


r  2tt?  2.~  2J2 

I d£  cos  — ^  exp  -  C  /2a  W  , 


(13a) 


a  NAEq  exp  -  2tt  a  f/X 


(13b) 


The  energy  transfer  for  a  bunch  is,  of  course,  proportional  to  the  energy 
transfer  per  particle  to  the  number  of  particles,  times  a  shape  factor.  The 
shape  factor  is  a  Gaussian  due  to  the  Gaussian  bunch  shape:  an  arbitrary 
localized  bunch  will  give  another  shape  factor,  dependent  on  the  same  para¬ 
meter  2ir  a 2  f/X.  In  any  case,  the  energy  transfer  does  not  decrease  too 
much  for 


a.  <  y  /X/TTf  .  (14a) 

The  allowable  bunch  width  aW  is  therefore 

aW  <  .  (14b) 

This  result  is  very  reasonable:  a  bunch  width  less  than  1/6  of 
the  free  space  wavelength  does  not  affect  the  energy  transfer.  Similar 
results  are  obtained  in  the  bunched  electron  accelerator  PHERMEX. 
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NUMERICAL  EVALUATIONS  AND  CONCLUSIONS 


The  energy  transfer  from  the  exact  cavity  fields  as  computed 
numerically  deviated  considerably  from  the  value  given  by  the  Hermite- 
Gaussian  approximation.  An  exact  formula  for  the  energy  transfer  is 
obtained  as  follows.  Consider  the  differential  equation 

[as-  h  *  VhV  -  r^-]  si„(h,n)  *  0  • 

multiply  this  by 

T  5  ell”//TV~  . 
and  integrate  from  n  =  -1  to  n  *  1 
Integration  by  parts  gives 

[ra-n2)sln  -  mV)  sj  = 

i 

jSln  [is  3n  *  Xln  -  hV  -  ITHt]  T  ^  • 

-1 

The  bracketed  operator  working  on  T  simplifies  considerably; 

£  Jt  -  (1-n2)’1*  C^ln-h2)  e  lhn  . 


s-lo 


A 


(15) 


(16) 


(17) 


(18) 


Therefore,  the  energy  transfer  integral  is 


]Vh^ 

J  nTTv 


dn  = 


:i-n2) (T'S-TS' 

X.  -h2 
In 


The  term  in  square  brackets  should  be  evaluated  in  the  limit  |n| 

In  this  limit  the  function  S,  is 

In 


sln(h,n)  =  Nln(i-n2)  k3h/2 


where  k^  for  large  h  is  found  as  (Slepian  196S) 


-h  n/2  (3n-l)/2 


k3  =  e  h  2 


♦  _ 

\T  ■ 


and  the  normalization  factor  N,  is  such  that  the  maximum  of  S,  is 

m  m 

approximately  unity: 


Nln  =  A/Cn-1)  !  tt/2 


The  square  brackets  in  (19) ,  exclusive  of  the  sinusoidal  term,  then  become 

simply  N  k,  h  .  The  energy  transfer  integral  becomes,  approximately, 
n  <5 


/  (n+s)/2  3n/2  -h  sin  .  ...  . 

U  *  ir\/(n-l)l  h  2  e  cos  n  /(-  ln'n  J 


(the  sin- function  should  be  chosen  for  n  even,  cos-  for  n  odd). 


This  result  should  be  used  with  caution,  because  the  estimate  for 

the  maximum  value  of  S,  contains  an  additional  constant  of  order  unity, 

in 
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Moreover,  the  above  formulas  use  only  the  first  terms  in  the  asymptotic 
expansions.  Figures  1-3  show  some  numerical  results. 


The  frequency  h  as  determined  from  the  WKB-approximation  to  the 
radial  spheroidal  wave  function  is  given  in  Figure  1:  In  Figure  1(a)  the 
axial  mode  number  n  =  1  ,  and  the  radial  mode  number  parametrizing  the 

curves  varies  from  Z  =  1  to  Z  =  21  .  The  data  are  in  general  agreement 

with  the  analytical  estimate  h  ^  (5^-1)^  (jj+...)  •  Figure  1(b)  given  the 

frequency  but  for  n  =  4  ;  The  frequency  is  higher  as  expected.  The 

frequencies  are  in  good  agreement  with  those  found  by  using  the  exact  radial 
spheroidal  wavefunction 

Figure  2  gives  the  parallel,  and  Figure  3  the  perpendicular  energy 
transfer.  Equation  (6),  for  l^n<5  as  function  of  the  cavity  position 
€0  and  parametrized  by  the  radial  mode  number  Z  .  Here  the  energy 
transfer  is  evaluated  with  Hermite-Gaussian  approximation:  Figure  2(a) 
for  n  *  1  shows  an  energy  transfer  that  increases  linearly  with  ; 

the  transfer  is  very  small  for  large  Z  ,  and  appreciable  for  Z  *  1  only. 


For  n  =  4  ,  as  in  Figure  2(b),  the  story  is  different:  Now  the 
energy  transfer  is  good  for  Z  =  1  and  EJ  — or  so,  and  for  larger  Z 
the  transfer  increases,  again  linearly  with  *  to  reasonable  values. 


Figure  2(c)  compares  these  data,  obtained  from  the  Hermite- 
Gaussian  approximation,  with  the  exact  formula  (19)  evaluated  numerically. 

Only  data  for  a  confocal  cavity  =  VT”  are  available  at  this  moment: 
the  figure  shows  AE|/eEpX  as  function  of  axial  mode  number  n  ,  parametrized 
by  Z  .  The'  surprise  here  is  that  the  agreement  is  not  uniformly  good: 
what  is  worse,  the  values  from  the  Hermite-Gaussian  approximation  are  too 
high,  by  a  large  factor. 


Probable  causes  for  this  discrepancy  are  under  investigation: 
one  possibility  is  indicated  by  the  exponential  factor  exp(-h)  in 
Equation  (23),  versus  exp(-h/2)  in  Equation  (6),  for  h  =  10  a  factor  of 
order  100  .  The  difference  in  exponential  decay  with  h  reflects  the 
unfortunate  fact  that  the  Hermite-Gaussian  approximation  is  not  good  in  its 
"tails";  the  axial  spheroidal  function  falls  off  more  rapidly  as  h  in¬ 
creases. 
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Figure  la 

The  frequency  h  =  ^  as  function  of  wall 
position  4q  for  axial  mode  number  n  =  1. 
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APPENDIX  V 

MISCELLANEOUS  ON  SPHEROIDAL  FUNCTIONS 


The  properties  of  the  spheroidal  functions  are  less  known  that 
those  of  the  more  common  special  functions.  In  this  section  various  useful 
formulas  are  given,  asymptotic  expansions  are  verified,  etc. 

Particularly  useful  sources  are  Meixner  and  Schafke  (Reference  1), 
and  Flammer  (Reference  2).  Numerical  data,  and  the  computer  code  used  by 
us,  were  obtained  by  courtesy  of  Dr.  A.  L.  Van  Buren,  NRL  (References  3 


* 


Normalization 


There  are  various  normalizations  in  use.  The  one  adopted  by 
some  modem  writers  is  Meixner's: 

I  [s-  (h’n) 

This  is  equivalent  to 


dq 


2  (n+m) ! 

2n+l  (n-m) ! 


p”  o’ 


where  P 


is  a  Legendre  polynomial.  The  set  S 


mn 


f(x)  =  2 

‘  n=m 


2n+l 


(n-m) 

(n+m) 


■:/« 


f(t)  S, 


mn 


is  complete: 
(h,t)  Smn(h,x) 


In  this  report  we  do  not  use  any  of  the  standard  normalization: 
instead,  the  function  that  gives  the  electric  field  is  normalized  to  unity 
at  its  maximum. 


APPROXIMATION  FOR  RADIAL  FUNCTION 


The  approximation  for  the  radial  function  R^Ch.Q  is 
(Reference  1,  Section  3.91) 

J.  <s)  •  (1) 

where 

s  =  ^  arc  cos  5_1  -  J  S2-l  ♦  •••  •  (2) 

Near  the  axis  5-1 


Here  q  =  2n-2m  +  l 

The  following  figure  gives  a  comparison  between  (52-l)* 
and  its  approximation.  For  frequencies  h  of  interest  h>10  , 

the  argument  is  extremely  close.  The  figure  shows  a  worst  case,  h  *  10 
and  n  *  5  .  The  "exact"  curve  is  the  left  side  of  Equation  (1),  the 

curve  marked  "appr"  is  the  right  side,  with  s  from  Equation  (2).  The 
dashed  curve  is  the  approximation  (1),  but  neglecting  the  last  term  in  s 
Each  function  is  normalized  to  the  maximum. 
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Comparison  Between  Exact  and  Approximate  Radial 
Spheroidal  Wave  Function 


ASYMPTOTICS  FOR  AXIAL  FUNCTION 


The  asymptotic  expansion  for  large  frequency  h-**-00  for 

the  angular  function  S  (h,n)  is: 

mn 

(l-n2p  S^fh.n)  =  Hn_m  (n/h)  exp  -  hn2/2 
00 

=  rS»  ar  Hr  exP  "  hn2/z 

where  Jl  =  n-ra  , 

o 

and  the  coefficients  a  satisfy  a  five-term  recursion  system  given  e.g. 

r 

in  Reference  2,  Section  8.  The  formulas  are  fairly  unwieldy:  The  following 
figure  is  a  comparison  between  (1-n2)  ’  Sln(h,n)  and  Hnl(rv\/h)  exp  -  hn2/2 
for  two  values  h  =  10  and  20  ,  and  two  values  n  =  2  and  4 

The  approximation  is  satisfactory  for  h  *  20  and  both 

values  of  n  ;  although  there  are  clear  deviations  for  n  =  5 
For  lower  frequency  h  =  10  the  agreement  is  acceptable  for  n  =  2  , 

but  bad  for  n  =  5.  Slepian  (Ref.  5)  gives  a  more  complete  treatment  of 
the  asymptotics  in  spheroidal  functions  which,  however,  is  fairly  complicated. 
The  relation  above  should  be  reasonably  good  in  the  region  |n|  £  h  ^ ,  but 
no  in  the  Gaussian  "tails"  where  n/~h  »  1. 
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V 

Comparison  Between  Exact  and  Approximate  Axial 
Spheroidal  Nave  Function  for:  h  ■  10  n  ■  2 


WKB  Approximation  for  Axial  Spheroidal  Function 


Quite  often  the  WKB  approximation  is  very  reliable,  especially  for 
the  higher  mode  envisioned  in  the  transverse  RF  accelerator.  To  our  knowledge, 
however,  this  approxmation  has  not  been  derived  before  for  the  spheroidal 
functions .  Consider  the  function 


un(h,n)  =  V  i-nz  sln(h,n) 


it  satisfies  the  Schroedinger  equation 


u"  s  Q  (*,n)u 


where 


q  =  HV-.A 

y  l-n2 

A  is  the  eigenvalue  for  the  spheroidal  function,  and  h  is  the 
frequency  h  =  fu/2c 


The  WKB -approximation  in  the  region  between  the  turning  points 
rit  *  y/\~/h  is 

r  nt 

u  —  2  (-Q)'*  sin  JyT- Q~dt  +  tt/4  .  (3 

n 

The  integral  can  be  done  exactly  in  terms  of  incomplete  Elliptic 
integrals.  Define  the  argument 


f  [”1-n*/n*l 

J  l-n2 
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the  integral  becomes 


Setting  n 


sin  0 


9t 

A  =  J  J  l-rf2sin20  d8 
9  « 


=  VT 


E(0t,n-1] 


E(8 


(5) 


where  =  arc  sin  nt 

The  symbol  E  is  the  incomplete  Elliptic  integral,  defined  by 

E(0,k)  a  j yi-k2sin20  d9  (6) 

Notice  that  in  Equation  (5)  the  parameter  nt_1  is  greater 
than  unity.  Existing  tables  are  typically  for  parameter  less  than 
unity.  For  computations,  therefore,  it  may  be  necessary  to  use  the  formula 
(AS  17.4.16) 

E(u|n~2)  =  n'1  E(Un^2,nt2)  -  Cn~2-i)u 


where 


E(U|m) 


o 


and  x  *  sn  U  .  Alternatively,  the  integral  can  be  evaluated  numerically 
as  it  is  needed. 


The  .approximation  for  the  spheroidal  function  Sln  becomes, 
between  the  turning  points 


5ln(h,n)  a  X"*  [(i-n2)(i-n2/nJ)]_J*  sin(vTA(nW4) 


(7) 
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while  the  approximation  for  the  derivative  of 


/7v~ sln  is 


_d 

dn 


'  -  „  IW/n2' 

yrVSi^b.n^ 


cos(vT-A+tt/4) 


(8) 


In  the  intermediate  region  about  the  turning  points,  where  the 

<v  / 

"potential"  is  linear  with  slope  the  approximation  is  in  terms 

of  Airy  functions:  with 


slB(h,n)  =  (l-n2)"5*  2 


,-7 


12 


(1-X/h2) 


*/i 


(9) 


In  the  region  outside  the  turning  points  •  nt  the  appropriate 
formula  is  similar  to  (3) : 


u 


Q  exp 


The  integral  now  becomes 


n 

/  /o”  dn 

nt 


(10) 


(11) 


it  can  also  be  written  in  terms  of  incomplete  elliptic  integrals,  but  that 
will  not  be  done  here.  Instead,  direct  numerical  integration  will  be  used 
when  necessary. 

Even  without  exact  evaluation  of  (11)  some  observations  can  be 
made  about  the  decay  rate  of  S^n  for  n  for  outside  the  turning 
point,  nt  ,  provided  nt  «  1  •  This  is  the  case  in  which  the 
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Gaussian-Hermite  approximation  applies  best.  \  comparison  between  the  two 
approximations  should  be  revealing,  especially  because  the  functions  occur 
in  the  energy  transfer  efficiency.  Slepian  (Ref.  5)  has  a  third  approximation, 
different  from  the  Gaussian-Hermite,  in  the  region  )nj  i  this  approxi¬ 
mation  should  agree  closely  to  the  WKB  expressions. 


J 
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APPENDIX  VI 


MISCELLANEOUS  ON  DEVICE  PHYSICS 


Current  Limitation  Revisited 


Maxwell's  equations  allow  a  simple  estimate  for  the  loading  of 
the  accelerator  cavity  by  the  beam  current.  Consider  the  equation 

V*?=  H  =  (1) 

c 

The  free  cavity  modes  calculated  earlier  in  terms  of  prolate  spherical 

functions  follow  from  (1)  by  setting  the  term  /tQj  to  zero.  If  fx  j  were 

-2 

everywhere  smaller  than  c  d E/  d t,  or  to  its  equivalent  VxB,  the  presence 
of  the  beam  would  scarcely  influence  the  wave  fields.  Let  us  compare  these 
terms  in  an  approximate  fashion. 

The  electric  field  vector  is  in  the  same  direction  as  the  current 
density.  It  is  clearly  negligible  if 

If  ^  (2) 

c 

Numerically,  this  means 

j  «  104  A/cm2  (3) 
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lOGhz. 


g 

for  an  accelerator  field  E  =:  10  V/m  and  frequency  v  =  ~ 

The  constant  current  density  can  be  spread  out  about  the  axis  over  an  area 
that  depends  on  on  the  acceptable  perpendicular  momentum  transfer,  or 
emittance.  Typically,  the  beam  electron  can  be  a  distance  X/20  from  the 
axis.  Therefore,  the  current  I  is  limited  by 


1  « 


E 

60/xov 


2 

10  A 


(3) 


for  the  given  parameters.  This  current  is  rather  small,  however,  for  a 

9  4 

field  E  —  10  V/m  and  v  =  lGhz  the  current  limit  already  10  A,  and 
currents  =  10^A  should  be  attainable  with  no  problems. 

Estimates  of  this  type  have  been  done  previously,  but  based  on 
energy  considerations,  with  the  same  results. 
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APPENDIX  VII 


DIFFRACTION  EFFECT  ON  MODE  STRUCTURE 


An  integral  equation  for  the  spheroidal  wave  function  is 
(Flammer,  5.3.12) 

1, 

ihCns 

•ln;h.n)  Rln(-  -  -  ■  i-  ■ 

-1 


x  J 


1-n  f 

)  =  i  J  ds 

-1 

1(htpT1flvY 


1-3  |Sln(l.,3)  • 


(1) 


Diffraction  due  to  open  cavity  ends  can  be  calculated  approximately  in  such 
an  integral  relation  by  changing  the  integration  limits  from  (-1,1)  to 
(-no.no)  ,  where  nQ  defines  how  large  the  cavity  hole  is. 


Equation  (1)  for  S^n  can  be  simplified  considerably  by  taking 
a  partial  derivative  3/3£  ,  and  setting  £  ■»  1  .  This  yields  first 

3Rln  1-n  [  ih?ns  /  £  / - 2  I - 7 

2SIn^T3i  Jdse  ^hnS  J1  *  ^  h^VVl-s2  Jx J 


sln(h.s)  ; 


* 


After  multiplying  by  -1  the  limit  £  ■»  1  can  be  taken,  with  the 
result 


lim 

5—1 


i~ 5 —  3R, 


In 


1 

1-n  f  ihns 
(h,n)  «■  i  Ids  e 

-1 


Note  that  the  limit  on  the  left  is  nonzero. 


This  formula  can  be  written  in  different  ways.  An  extra  factor 


^'l-ri  gives 


i  An  Un(h,n)  =  (l-n  )  Ids  e  un(h,s)  , 


where 


unCh,n)  =  sln(h. 


n) 


n  K 


lira  X2  .  dRln  2 

— ~1  3£  h 

m  “ 


Comparing  Equation  (4)  to  the  well-known  analogue  for  the  zero 
order  spheroidal  function 


2i  R  (hfl)  S  (h,n)  *  /  e  S  (h,s)ds 

on1  on  I  on 


we  note  an  extra  factor  (l-n  )  •  In  addition  the  eigenfunction  is  not  a 

_ • 1 1 :  1  _ i  c _ _ 1 .  1 r* „  . f  .  _2  _ .  . 


purely  spheroidal  function,  but  has  an  extra  factor  \  l-n  ,  to  make  it 
the  potential  for  the  normal  electric  field 


E^  « 


3Un(h,n) 


The  interpretation  of  Equation  (4)  is  then  that  the  normal  electric  field 

2 

pattern  on  the  cavity  wall  is  reflected  with  the  propagator  (l-n  )  exp  ihn* 


Other  forms  of  Equation  (4)  are: 


n-1 

i  An  SlnCh,n) 


=  ^1-T|2  jds  e  ^  ^1-s2  Sln(h,s)  ds  , 


which  would  propagate  the  ^component  of  the  magnetic  field,  or 


1  An  Vh"n)  - 


'iy4-  /ds  e  W(l-s2)(e^-s2)  Vn(h,s) 
£  -n  J  ” 

0  -1 


•  (10) 


where 


Vn(h,n)  = 


Sln^ 

V_2  2 

Vn 


This  operator  propagates  the  parallel  electric  field 


E  «  V  (h.n)  • 

n  n 


Natice  that  Equation  (10)  seems  to  contain  the  cavity  wall  position 

£0  explicitly  (not  only  through  the  frequency  h  ) .  However,  the  Equations 

(4)  and  (9)  show  that  this  is  not  the  case:  the  explicit  £q  -  dependence 

is  cancelled  by  the  extra  £  hidden  in  V  (h,n)  of  Equation  (11).  Since 

o  n 

ail  three  integral  forms  are  equivalent  any  of  these  can  be  taken  as  a  basis 
for  perturbation.  On  additional  form  can  be  obtained  by  setting  £  ■  1  in 
in  (10)  and  (11): 


where 


L  Jr  ihns  - 

An  vn(h»n)  -  /  ds  e  (lV)  V(h, 


s) 


(13a) 


(13b) 


Diffraction  effects  on  the  cavity  fields  are  calculated  most 
conveniently  with  Equation  (9),  by  changing  the  integration  limits  from 
(-1,1)  to  (-n0,n0)  •  The  wave  structure  in  the  open  cavity  X(h,n)  is 

now  the  eigenfunction  of 

% 

n-i  / - Y  f  ihns  / - 2 

i  B  X(h,n)  =  Jl-n  /  ds  e  Jl-s*  X(h,s)  ,  I 


where  B  =  Bn(nQ)  is  the  eigenvalue.  >btice  that  B  is  purely  real  in  the 

case  n  — "*1  ,  B  =  A  ,  and  that  in  this  limit 

o  n 


X(h,n)  =  xn(h.n;n0)-^  sln(h,n)  ; 


Therefore  X  has  the  same  number  of  zeroes  and  symmetry  properties  as  S^n  ; 
even  for  n  -  1  even,  and  vice  versa. 


For  small  An  =  1-n  the  integration  in  (14)  can  be  approximated 
by  subtracting  the  average  value  over  the  tail  ends  (-l,-no)  and  (nQ»l) 
from  the  standard  integration  interval  (-1,1)  in  Equation  (9).  The 


result  is 


n-l  f  ~  Y  ihns  I  Y 

i  B  X(h,n)  =  Vl-n  Ids  e  e/l-s*  X(h,s) 


n-l  3/2  | - r-  cos,.  . 

-  2i  (An)  X(h,na)Wl-n  sininnV  .  (16) 


Here  cos(hnn  )  applies  when  n-l  is  odd;  n  is  a  average  between 

&  3L 

1  and  n0  »  na  *  (l-»n)/2  ,  and  An  -  i-nQ  • 
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The  perturbation  term  is  proportional  to  (An)  ,  due  to  the 
approximation  of  the  square  root  under  the  integral.  Also,  the  perturbation 
has  the  same  symmetry  properties  as  the  left  side  of  (16),  and  the  same 
factor  in  *  .  Therefore  this  first  order  perturbation  will  not  produce 

an  imaginary  part  in  the  eigenvalue  B  . 

Following  straightforward  perturbation  theory  the  unknown 
function  X(h,n)  can  be  expanded  in  a  series  of  spheroidals: 


X(i.,n)  =  Jj  an  slnCh.„) 


(17) 


where  the  coefficients  are  given  by 

■  si Sit  sln(h.n)dn  •  (18) 

-1 

The  constant  (n+!s)/n(n+l)  is  appropriate  to  the  Meixner  normalization 
of  the  spheroidal  functions 


(n-wi) ! 

(n  +^)  (n-ra)  ! 


(19) 


Equation  (16)  becomes,  after  interchanging  integration  and 
summation,  and  using  the  defining  integral  (9) : 


2 
•  n 


(anB  - 


a  A 
n  n 


Cn> 


sinCh.n) 


0 


(20) 


where 


c 

n 


, An1*3/2  Y  i-K  n  2k  +  l 

(An)  x(h,na)  k£k+1j 


)  P"  sin 


hnna  Sln(h,n^ 


(21) 


Evaluation  of  X(h,n)  demands  another  series,  but  to  lowest  order 


x(h,na)  •  sln(h.n4) 


(22) 


Since  the  S^'s  are  orthogonal  there  follows 

an  *  A~-B  • 
n 


(23) 


The  perturbed  eigenvalue  B  is  obtained  from  approximating  X(h,n)  =  X^h.njr^) 
in  Equation  (18)  by  Sin(hjT0  .  Then  aN  *  1  311(1 


For  practical  calculations  the  coefficients  must  be  evaluated 
but  this  is  quite  tedious.  An  alternative  approach  is  to  solve  Equation  (9) 
directly  by  computer.  A  suitable  algorithm  is  the  interation  of  Equation  (4) . 


Having  found  the  fields  as  function  of  n  inside  the  open  cavity, 
from  a  perturbation  series  or  by  iteration  of  Equation  (14),  it  is  straight¬ 
forward  to  calculate  the  5  ~  dependent  factor  in  the  fields,  Rln(h,£) 

The  derivative  of  Equation  (1)  with  respect  to  n  ,  an  extra  multiplication 
with  '  ,  and  the  limit  n^l  yields 
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in  complete  analogy  to  Equations  (3)  and  (4) .  However,  now  the  eigenvalue  is 


C 

n 


/,  2  3S 

V1"1  5? 

lim  n-*-lJ 


(26) 


For  an  open  cavity  the  ^-dependent  factor  Yn(h,£)  is  now  deter¬ 
mined  by  changing  the  integration  limits  in  (25)  from  (-1,1)  to  (-r)  ,nQ). 

It  is  not  immediately  obvious  that  this  is  the  correct  way  to 
find  the  £ -dependence,  but  this  algorithm  is  certainly  reasonable:  it 
makes  the  field  Ez  continuous  across  the  focal  point  0=1,  £=1.  To  verify 
the  continuity  one  should  remember  that  along  the  axis  has  two  distinct 
forms. 


E  =  E  between  foci 

z  o 


E 


z 


outside  the  foci 


(27) 


The  explicit  expressions  for  E^  and  evaluated  at  n  =  £  s  1  then  show 
the  continuity  at  this  point. 


The  field,  parallel  to  the  cavity  axis,  Ez  *  E^,  evaluated  on 
the  axis  £  =  1  between  the  foci  is: 


i 


h 


ASln(h,ti) 

P 


• 

3 

5C  \ 

A2-1 

. 

S=1 


(28) 
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Here  A,  which  depends  on  the  normalization  of  the  various  spheroidals  is 
a  field  strength  parameter,  whereto  we  shall  return. 


Inserting  Equation  (9),  with  the  explicit  value  (3b)  for  its 


coefficient  An  gives 


iA 


tyS-i.n)  =  -  ~ 


R 


In 


\T*'h 


In 

C-l 


Tn(h,n) 


(29) 


where 


Tn(h'n)  =  /  ds  elhsn^/l-s2  Sln(h,s) 


(30) 


On  the  other  hand,  E  =  E-  outside  the  foci.  On  the  axis  (n=l,5) 

Z  s 


the  field  is  explicitly 


Ain  Rm(h'^ 


UC  h 


3 

3n 


Sln(h,n) 


(31) 


mi 


Use  of  Equations  (25)  and  (26)  for  the  function  Rln  gives 


EE(?.n-i)  ■  -  ~ 


fnV*3 


2  e 

n  s 


in 


2  3_fln 
3n 


T„(h^) 


(32) 


ml 


8-8 


The  field  Ez  is  continuous  across  the  focal  point  £  =  1  =  n 
provided  that  the  constants  in  the  square  brackets  are  equal.  This  is 
easily  seen  directly  from  the  integral  relations  (9)  and  (25).  Namely, 
the  square  brackets  do  not  depend  on  the  normalization  constants  An  and 
Cn,  but  are  simply  functionals  of  Tn(h,n)  and  Tn(h,£):  and  it  is  clear 
from  definition  (30)  that 

Tn(h,n  -  1)  =  Tn(h,5  -  1)  ,  (33) 

Moreover,  this  equality  is  true  for  all  derivatives. 

Numerical  results  are  shown  in  Figure  1  and  2.  Figure  1  is  the 

field  along  the  z-axis  for  the  case  h  =  6,  n  =  2,  n  between  0  and  1,  and 

£  from  1  to  3.  The  parameter  is  the  size  of  the  cavity  opening,  parametrized 

by  q  as  indicated.  The  fields  for  the  closed  cavity  are  concentrated  most 
o 

between  the  foci,  but  as  the  cavity  opens  up  the  field  pattern  widens,  and 
RF  leakage  increases. 

The  coupling  from  RF  fields  to  particles  generally  deteriorates 
from  opening  the  cavity.  Figure  2  shows  the  coupling  for  the  case  h  =  6,  n  =  2 
as  function  of  cavity  hole  size  nQ,  for  the  lowest  three  radial  modes 
i  =  1,  2  and  3.  The  coupling  is  not  affected  too  much  for  the  lowest  mode 
up  to  nQ  a  0.8.  The  higher  modes  are  more  sensitive,  to  the  point  of 
reversing  the  sign  of  the  coupling  constant. 

The  conclusion  from  the  diffraction  calculations  is  that  the 
cavity  can  be  safely  opened  somewhat,  e.g.  up  to  nQ  a  0.9,  without  too  much 
decrease  in  the  coupling.  However,  there  is  no  increase  in  coupling  efficiency 
by  increasing  leakage  from  the  cavity. 
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Figure  1.  The  z-component  of  the  electric  field  at 
frequency  h  ■  6,  and  axial  mode  number 
n  -  2,  along  the  cavity  axis  for  an  open 
cavity. 

The  parameter  n0 

nQ  ■  1  -  0.05  m,  m  ■  0  -  8  (except  nQ  *  0.95). 
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Figure  2.  The  coupling  from  RF  fields  to 

relativistic  particles  as  function 
of  the  cavity  opening  for  three 
radial  modes  and  h  »  6,  n  *  2. 
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